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ABSTRACT 

Context. Sh2-104 is a Galactic Hn region with a bubble morphology, detected at optical and radio wavelengths. It is considered the first observa- 
tional confirmation of the collect-and-collapse model of triggered star-formation. 

Aims. We aim to analyze the dust and gas properties of the Sh2-104 region to better constrain its effect on local future generations of stars. In 
addition, we investigate the relationship between the dust emissivity index {S and the dust temperature, T dust . 

Methods. Using Herschel PACS and SPIRE images at 100, 160, 250, 350 and 500 //m we determine T dlL „ and /? throughout Sh2-104, fitting the 
spectral energy distributions (SEDs) obtained from aperture photometry. With the SPIRE Fourier transform spectrometer (FTS) we obtained spec- 
tra at different positions in the Sh2-104 region. We detect ./-ladders of 12 CO and 13 CO, with which we derive the gas temperature and column 
density. We also detect proxies of ionizing flux as the [Nil] 3 _P, — 3 P and [Ci] 3 P2 — 3 Pi transitions. 

Results. We find an average value of fi -1.5 throughout Sh2-104, as well as a T dust difference between the photodissociation region (PDR, ~ 25 K) 
and the interior (~ 40 K) of the bubble. We recover the anti-correlation between /? and dust temperature reported numerous times in the literature. 
The relative isotopologue abundances of CO appear to be enhanced above the standard ISM values, but the obtained value is very preliminary and 
is still affected by large uncertainties. 
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1. Introduction 



Sharpless 104 (Sh2- 1 04. [Sharplessl[T959l) is an optically visible 
Galactic Hn region with a bubb l e morphology, excited by an 
06V star dCrampton etall Il978t iLahullal Il985h . It is located 
~ 4kpc from the Sun dDeharveng et al.Ll2003l) . with galactic co- 
ord inates 74.7620; + 60 (J2 000). 

iDeharveng et alJ (120031) proposed Sh2-104 as a strong 
candidate for massive trigge red star formation throug h the 
collect-and-collapse process dEhnegreen & Lada 1 19771) . The 
ionize d region is also visible at radio wavelengths (Fich, 
1993), and an ultracompact (UC) Hn region, coincidant with 
the IRAS 20 160+3636 source, lies at its eastern border 
(ICondonetal.lll998h . 

We present new submm images and spectra taken towards 
Sh2-104 with the Herschel Space Observatory dPilbratt et al.L 
2010). These observations allow us to map a wavelength range 
not easily accessed before, providing new insights into the dust 
and gas properties of Sh2-104. 



"Evolution of Interstellar Dust" of SPIRE dAbergel et al., 2010). 
andHOBYS of PACS ( Motte et alll20Toh . A 15'xl5' region was 
imaged with PACS at 100 and 160 /im (at resolutions of 10" and 
14"), and with SPIRE at 250, 350 and 500 //m (at resolutions of 
18", 25" and 36"). Spectra were taken with the SPIRE-FTS long 
and short wavelength receivers (SLW and SSW, respectively) 
at seven different positions with sparse sampling, covering the 
194-671 /mi range. The resolution at the receivers' central pix- 
els varies between 17 - 19" for SSW and 29-42" for SLW. The 
data were reduced with the HIPE software version 2.0 with the 
latest standard calibration dSwinvardet al.. 2010). 

Figure [TJ shows a color-composite image of Sh2-104 with 
PACS 100 /Km (blue), SPIRE 250 //m (green) and SPIRE 500 yum 
(red). Different regions of interest, addressed in the following 
sections, are superimposed. We can see that the interior of the 
bubble is brighter in the PACS band, showing the hotter temper- 
atures of the material in this region. On the other hand, outside 
the bubble the material is colder and emits stronger in the SPIRE 
bands. 



2. Observations 

The Herschel observations were taken on 2009 D ecember 17 
simult aneously with PACS dPoglitsch et all 1201 Oh and SPIRE 
(Griffin et al., 2010), as part of the guaranteed-time key-projects 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



3. Dust properties 

We assume the dust emission in Sh2-104 can be modeled by 
an (optically thin) gray-body and that the emissivity of the dust 
grains can be fitted with a power law 



S v oc QB v {T)k()\ — 



N. 



(1) 
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Fig. 1. Composite image of Sh2-104. The field is ~ 13' x 13' and 
shows PACS 100/im emission in blue, SPIRE 250 //m in green 
and SPIRE 500 pm in red. Outlined are the regions where aper- 
ture photometry was applied, the region used for background 
substraction is off the map. Blue and red circles mark the point- 
ings of the central pixels of SPIRE-FTS SSW and SLW, respec- 
tively. Their diameter represents their respective average resolu- 
tion of 19" and 35". 



where S v is the measured flux density at frequency v, Q is 
the observing beam solid angle, B v (T) is the Planck func- 
tion for temperature T at frequency v, kq (v/vq)^ is the dust 
opacity, j3 is t he dust emi s sivity index, and N is the dust col- 
umn density (Hildebrand, 1983). The value of p is believed 
to range between ~ 1 and ~ 2, but is an open issue that is 
still discu ssed, as is its dependence on the dust t emperature 
(see e.g.,|O ssenko pf & Henninglll994T) . lDupac et ail d2003l) and 
IDesert et all J2008) foun d an inverse relation ship between/? and 
Tdust- On the other hand. lShettv et alJ t2009) suggested that this 
result arises from noise and the combination of multiple emis- 
sion components along the line of sight. 

To determine the dust temperature structure of Sh2-104, we 
performed aperture photometry measurements on selected areas 
in the field. The apertures are shown in Fig. Q] and sample the 
interior of the bubble and the photodissociation region (PDR), 
including the UCH n region associated with IRAS 20160+3636. 
We used a single aperture (not shown) to account for the back- 
ground emission. We fitted the PACS+SPIRE emission for all re- 
gions. These data represent the cold emission component, there- 
fore we fitted a single temperature. We did this first allowing p 
to vary, and later fixing p =1.5. The resulting Tdust and/3 values 
are shown in TableQ] their uncertainties are the formal 1 cr values 
from the fitting procedure. Figure|2]shows examples of the fitted 
spectral energy distributions (SEDs) for the regions UCHII and 
PDR 1. Apart from the clear difference in Tdust, it can be seen 
that for UCHII the PACS+SPIRE emission does not allow to 
sample the peak of the SED, which is reflected in a larger error in 
the fit. This is also seen in the four Interior # regions, suggesting 
the presence of a warmer emission component that contributes 
to the shorter- A emission. Observations in the mid-IR can pro- 
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Fig. 2. Example of the SEDs obtained for regions UCHII (filled 
red squares) and PDR 1 (open blue squares). A dashed line is the 
fit with p = 1 .5 and a solid line is the fit with p allowed to vary, 
with the resulting (TdusuP) values shown in the respective color. 
The temperature difference between the regions is apparent, as 
is the difference in the peak sampling of the SED for PDR 1 and 
UCHII. 



vide a constraint on this component, and a 2-temperatures fitting 
would be more appropriate to determine the Tdust and/? values. 

The temperatures throughout the PDR are between ~ 20 
and ~ 30 K, and the UCH n region is marginally warmer. The 
region outside the bubble (region Outer) is the coldest, while 
the regions mapping the interior of the bubble are hotter, with 
an average temperature of ~ 40 K. The average /3 obtained is 
~ 1.5, therefore the temperatures obtained from the fit with a 
fixed p - 1 .5 do not significantly differ from those obtained with 
a free/? (see Table[TJ. 

Figure [3] shows the distribution of the spectral indices p 
vs. dust temperatu re, a long with the rela tionships found by 
iDupac et all (120031) and IDesert etaf] d2008l) . Within the uncer- 
tainties, our results agree with both relationships. We can also 
identify two groups, one with temperatures between 35 and 47 K 
and an emissivity index between 1.0 and 1.3, and the other with 
T - 17 -29 K and/? = 1.5 - 1.8. These two groups show that on 



Table 1. Results from aperture photometry. 



Region 


Tdust (K) 

p 


P 

-ree 


Tdust (K) 

P= 1.5 


CS 


25 ±3 


1.5+0.3 


25.5 ± 0.9 


Interior 1 


43 ± 12 


1.2 + 0.4 


35.7 ±2.1 


Interior 2 


39 ±9 


1.3 + 0.3 


33.5 ± 1.7 


Interior 3 


47 ± 15 


1.0 + 0.4 


34.2 ± 1.9 


Interior 4 


35 ±8 


1.3 + 0.4 


30.4 ±1.4 


Outer 


17 ±2 


1.7 + 0.4 


17.3 ± 0.4 


PDR1 


20 ±2 


1.7 + 0.4 


21.7 ±0.6 


PDR 2 


22 ±3 


1.7 + 0.3 


24.1 ±0.8 


PDR 3 


23 ±3 


1.7 + 0.4 


24.6 ± 0.8 


PDR 4 


22 ±3 


1.5+0.3 


22.3 ± 0.7 


PDR 5 


22 ±3 


1.5+0.3 


22.9 ± 0.7 


PDR 6 


24 ±3 


1.8 + 0.3 


26.0 ± 0.9 


PDR 7 


25 ±3 


1.7 + 0.3 


26.9 ± 1.0 


PDR 8 


26 ±4 


1.7 + 0.4 


28.6 ±1.2 


UCHII 


29 ±5 


1.8 + 0.3 


32.7 ± 1.5 
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Fig. 3. Distribution of the emissivity index /3 vs. dust temper- 
at ure. The dotted and das hed lines are the re lationships found 
bv lDupac etal] d2003h and lDesert etal] d2008l) . respectively. We 
cannot distinguish between the two relationships, but lower /3 
values correspond to higher temperatures (see text). 
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Fig. 4. SPIRE-FTS spectrum corresponding to the central pix- 
els of SLW (pixel C3, red) and SSW (pixel D4, blue), towards 
the UCHii region associated with IRAS 20160+3636 (pointing 
E). The lines detected are labeled. The jump between continuum 
levels is a calibration effect. 



average higher/? values are preferentially associated with colder 
material. 

Although anti-correlation between B and Tdust is reported 
in the literature (e .g., iDupac et alj|2003l: lYang & Phillipsll2007t 
iDesert et al1l20 08). it is yet not clear which physical processes 
are behind it. Other authors studying this relationship examine 
the emission from different regions scattered in the sky. In con- 
trast, we are finding this anti-correlation in the analysis of one 
contiguous complex object at a specific location in the sky. A 
f5 - Tdust relationship ma y indicate a change of the dust proper- 
ties ([Stepnik et al., 2003). Grain fluffiness in particular increases 
the e missivity index while keeping a relatively lo w temperature 
(e.g.. lStognienko et al]ll995UFogel & Leundll998l) . 

Fluffy grains result from grain coagulation and growth. The 
grain coagulation timescale and feasibility depend on factors 
such as the existence of ice mantles, grain size, and relative grain 
velocities. 

We are finding the highest /3 values and lowest Tdust values 
toward the PDR of Sh2-104, which would imply then that the 
fluffiest and largest grains are located in the PDR. The question 
remains whether the PDR dust coagulated after the creation of 
the H H region, or if the birth of the H n region has destroyed the 
already coagulated dust located in the ionized cavity. 

4. Gas properties 

The seven pointings of the central pixels observed with SPIRE- 
FTS are marked with red (SLW) and blue (SSW) circles 
in Fig. [T] In total we detected the transitions described in 
Table [2] The richest spectrum (Fig. |4]i was obtained towards 
pointing E, which targets the UCHn region associated with 
IRAS 20160+3636. The most prominent features are the CO 
/-ladder and the [Nn] 3 P; - 3 Pq transition, an importa nt ionized- 
gas c oolant and a proxy for the Ha flux (see e.g., lOberst et al.l 
120061) . 

In the simple hypothesis of optically thin emission, we plot- 
ted 12 CO and 13 CO exci tation diagrams, following the formu- 
lation of I Johansson et al.l (1198 4). An example is shown in Fig. 
[5] for pointing E. The slope of the linear fit for I2 CO results in 
a temperature of T n = 246 + 2K, and T n = 170 + 53 K for 



13 CO. The total column densities derived are TVJ 2 = 10 16 • 36±001 
, respectively. The fits in Fig.[5]do not 



and/V| 3 = 10 15 - 4±01 cirT 2 



include all the measurements. The downturn seen in the lower 
levels of 12 CO is interpreted as the optically thick/thin regime 
turnover, and is most likely a real physical feature and not an 
instrumental or calibration effect, because it is only seen for that 
species and not, for example, for 13 CO. If the lines were emitted 
from an optically thick medium, their intensity would be under- 
estimated, thus their respective column densities would also be a 
lower limit. Therefore, the optically thin assumption would hold 
for 12 CO only for transitions higher than / = 9 — > 8, and those 
are the points included in the fit. 

For 13 CO on the other hand only the lines with wavelengths 
in the SPIRE-FTS SLW range are detected. These correspond to 
the 7 = 5->4 to 7 = 9->8 transitions. The 13 CO lines in the 
SPIRE-FTS SSW wavelenght range (J = 10 -» 9 to J = 14 -> 
13) are detected as upper limits, because the line positions are 
found displaced from their expected positions, indicating that we 
are probably seeing some "outlier" noise features rather than the 
lines themselves. Therefore, we fit and show in Fig.[5]only the 
five 13 CO transitions detected with SPIRE-FTS SLW. Following 
the reasoning of the previous paragraph, the 13 CO transitions are 
most likely optically thin. 

The two distinct gas temperatures obtained with 12 CO and 
I3 CO suggest two different gas components or a stratification 
of the emitting region. The colder component is traced by the 

Table 2. Lines detected with SPIRE-FTS in the different point- 
ings. 



Transition 


Rest A fjum) 


E,„, (K) 


Transition 


Rest A (fim) 


E„„ (K) 


CO(4 - 3) 


650.3 


55.4 


1J CO(5 - 4) 


544.2 


79.4 


CO(5 - 4) 


520.2 


83.0 


13 CO(6 - 5) 


453.5 


111.1 


CO(6 - 5) 


433.6 


116.3 


13 CO(7 - 6) 


388.7 


148.2 


CO(7 - 6) 


371.7 


155.0 


13 CO(8-7) 


340.2 


190.5 


CO(8 - 7) 


325.2 


199.3 


13 CO(9 - 8) 


302.4 


238.1 


CO(9 - 8) 


289.1 


249.1 


[Ci] 3 />, - 3 P 


609.1 


23.6 


CO(10 - 9) 


260.2 


304.4 


[Ci] 3 P 2 - 3 Pj 


370.4 


47.3 


CO(ll- 10) 


236.6 


365.3 


CH + (1 -0) 


359.0 


40.1 


CO(12- 11) 


216.9 


431.7 


[Nil] 3 P, - 3 P 


205.2 


70.2 


CO(13- 12) 


200.3 


503.6 
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Fig. 5. 12 CO (red squares) and 13 CO (blue diamonds) rotational 
diagram for pointing E towards the UCHn region. The solid 
lines are the best linear fit. Points included in the fit are shown 
as filled, points excluded as open. The turnover into the optically 
thick regime is noted in the lower-level transitions of 12 CO. 



optically-thin 13 CO transitions at energy levels for which 12 CO 
is optically thick, while the hotter component is traced by the 
more energetic, optic ally- thin 12 CO transitions. Therefore it is 
likely that the colder gas is located at grater depths in the PDR 
than the hotter gas. 

Assuming similar emitting volumes and beam filling fac- 
tors as well as optically thin emission, our CO and 13 CO col- 
umn density values imply an abundance ratio [ 12 CO]/[ 13 CO], 
which is several ti mes lower t han the reported elemental value 
of 12 C/ 13 C ~ 69 dWilsonL[T999]) . This would imply an enhance- 
ment of the 13 CO isotopologue abundance. 

However, several factors can contribute to the low 12 C/ 13 C 
abundance ratio we find. Perhaps the most important one would 
be the assumption of optically thin emission for 12 CO. We used 
the high-energy transitions to derive its column density, and al- 
though in a first analysi s they appear to be o ptically thin, it might 
not be the case (see e.g jHabart et alj2010t) . To address this issue 
we will present and analyse PDR models of the Sh2-104 region 
in a forthcoming paper. 



5. Summary 

With Herschel PACS and SPIRE data we have analysed dust 
and gas properties of the bubble-shaped Hn region Sh2-104. 
Aperture photometry of PACS+SPIRE images allowed us to 
derive the dust emissivity index j3 and the dust temperature 
throughout Sh2-104. We found two different groups, one at 
colder temperatures and higher j3, and the other at warmer tem- 
peratures and lower j3. We recover the inverse ft - Tdust rela- 
tionship reported in the literature, but the uncertainty in our 
fit prevents us from finding a precise depenc e. As an e xam- 
ple, the diff erent relationships found by Dupac et al. (2003) and 
iDesert et a l. (2008) are both compatible with our results. PACS 
and SPIRE have allowed us to constrain the "cold" side of the 
SED. We have estimated both/? and Tdust simultaneously, but the 
uncertainties remain relatively high. With the data at hand we 
can only fit a single temperature component, which disregards 
the contribution of warmer components to the PACS shorter 
wavelengths. The existance of a/3 - Td ust anti-correlation could 



be due to differences of the dust grain properties between the 
PDR and the H n cavity. 

SPIRE-FTS spectra at different pointings throughout Sh2- 
104 have unveiled the CO chemistry of the region in more detail. 
We detect the 12 CO and 13 CO /-ladders up to the J = 13 -> 12 
and ./ = 9 — > 8 transitions respectively, revealing the warm 
gas component in the region. Rotational diagrams towards the 
UCH ii region in the PDR of Sh2-104 show that the i3 CO emis- 
sion is optically thin and also the 12 CO for transitions above the 
7 = 8 level. The emission shows two different gas components, 
a colder one with a temperature of ~ 170 K and a hotter one at a 
temperature of ~ 250 K. The CO column densities derived would 
suggest an enhancement of the 13 CO isotopologue abundance ra- 
tio with respect to the elemental value, but the uncertainties of 
the different assumptions are still too large to confirm that re- 
sult. In a follow-up paper we will show models of the PDR of 
Sh2-104, which will provide better constraints on the gas tem- 
perature, density and column density structure. 
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